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Abstract 

We study an S-I type epidemic model in an age-structured population, 
with mortality due to the disease. A threshold quantity is found that 
controls the stability of the disease-free equilibrium and guarantees the 
existence of an endemic equilibrium. We obtain conditions on the age- 
dependence of the susceptibility to infection that imply the uniqueness 
of the endemic equilibrium. An example with two endemic equilibria is 
shown. Finally, we analyse numerically how the stability of the endemic 
equilibrium is affected by extra-mortality and by the possible periodicities 
induced by the demographic age-structure. 

1 Introduction 

In this paper we study the existence and uniqueness or multiplicity of positive 
steady states for an S-I type epidemic through an age-structured population. 
The population growth is subject to intra-specific mechanisms that, in the ab- 
sence of infection, lead to a steady state which can be destabilised by Hopf 
bifurcation, for some values of the intrinsic reproduction number Rq. The in- 
fection is regulated by an age-dependent contact mechanism and interferes with 
the demographic process introducing additional mortality due to the disease. 

A model with disease-induced mortality was for the first time introduced in 
PQ and the formulation in an age-structured context is due to May in [IB]. In 
[2] Andreasen analysed two different cases: an age-structured epidemiological 
model and a model where infection has a constant duration. For the former 
he assumed that only the epidemiological interactions are density-dependent; 
he obtained conditions for the existence of one endemic state and studied its 
stability in the hypothesis that only fecundity is age dependent. 

Here we give conditions for existence and uniqueness of an endemic equi- 
librium. A situation where two endemic equilibria appear is shown. Multiple 
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solutions are specific of this model, in the sense that they do not occur in mod- 
els without age structure or in the same model without the extra-mortality. 
Models with multiple steady states had already been proposed in the literature, 
the current model exhibits a different mechanism through which they arise. In 
particular, in [TJ] the role of variable population size in a multigroup epidemic 
model is emphasised, showing that multiple endemic equilibria are possible. 
Concerning the analysis of epidemic models with multiple endemic equilibria in 
an SIS epidemic model without age structure, see also [TS] and the bibliography 
therein. 

In addition to the study of existence of equilibria, we also investigate their 
stability showing that the effect of disease-induced mortality may produce sta- 
bility changes related to Hopf 's bifurcation. Using a numerical method exposed 
in [3], we explore the behaviour of our model, within some critical parameter 
regions. 

The paper is structured as follows. Next section is devoted to present our 
model, while in Section [3] the well-posedness of the system of partial differen- 
tial equations is pointed out. In Section @] conditions for existence of endemic 
equilibria are found. Successively, in Section [5] uniqueness versus multiplicity of 
endemic states is discussed. 

In Section [5] the characteristic equation is computed, so that in Section[7]the 
stability analysis of the disease- free equilibrium is carried out and in Section [5] 
a stability change is described. 

Finally, in Section [9] stability is numerically analysed in two different situa- 
tions: where multiple endemic states occur and where periodic solutions arise 
through Hopf 's bifurcation. Our analysis reveals the importance of the param- 
eter of disease-induced mortality. In the last section conclusions are drawn. 



We consider a population that, in the absence of infection is described by N(a, t), 
the age density at time t > 0, where a £ [0, a'], < +oo being the maximum 
age an individual of the population may reach. The growth of the population 
is regulated by the following model of Gurtin-MacCamy type (see [10) 1 



In (|2.ip . /3(a) and fJ-(a) are the intrinsic vital rates, with /3 normalised to 



2 The model 



— + — +M(a)JV(a,t) =0, 



< 




(2.1) 



satisfy 




where 



jt(o) = e-ti^ 
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is the survival probability, i.e. the probability at birth of surviving to age a. 
Moreover, the parameter Rq is the demographic basic reproduction ratio and 
we suppose that 

Rt > 1 . 

The function $ is a Lipschitz-continuous non-increasing function describing 
density-dependence of births. We assume that 

$(0) = 1 , lim $>(x) = 

x — >+oo 

and that there exists xo £ M such that &(xo) < V-^o an< ^ ^KO ^ s decreasing on 
[0,a;o]- Finally, Q is the size 

r 1 ' 

Q(t)= / r(a)N(a,t)da, (2.2) 
Jo 

where r(a) is a weight kernel. We point out that, by the assumptions made, 
equation 

Rfo(Q d ) = 1 (2.3) 
has a unique solution Q* d , providing a non trivial equilibrium 

N * {a) = ( 2 .4) 

Jo r(a)n(a)da 

of the population. 

We recall (see [13]) that, provided that $ is differentiable at Q d , the stability 
of this steady state is related to the characteristic equation 

e- Xa p{a)Tr(a)da + R^(j)'(Q* d )N*{0) / e- Xa r{a)n(a)da , (2.5) 

Jo 

the analysis of which also provides information about possible bifurcation points. 
For the problem above we make the following standard hypotheses: 

/?(•) £ L°°(0, a 1 ), 0(a) >0 in [0,0*], 
M(-) 6 4, c (0,ot), ^(o) >0in [0,ot], 

r(-)6i"(0,a f ), r(a) >0 in [0,a f ]. 

Since no individual may live past age a', in order to have 

jt(o + ) = 0, 

we need to assume 

,.t 

/ mC 17 ) ^ cr = +°° ■ 

Jo 



3 



The population develops an S-I type epidemics so that, denoting respectively 
by S(a,t) and I(a,t) the age-specific densities of susceptible and infective indi- 
viduals at time t, the following system of equations describes the transmission 
dynamics of the disease: 



as dS_ 

dt da 

dl dl 

dt da 



---(\(a,t)+n(a))S(a,t), 
X(a, t) S(a, t) - {n(a) + a) I(a, t) , 



(2.6) 



S(0, t) = Rfo(Q(t)) / /3(a) (S(a, t) + I(a, t)) da , 
Jo 

1(0, t) =0, 
S(a,0) = S (a) >0, 
{I(a,0)=I (a) >0. 

Here the terms A(a, t) and a describe the mechanism of infection. In particular, 
a is the extra-mortality due to the infection and A(a, t) is the force of infection. 
Concerning Q(t), in (|2.6|) we set 



Qit) = / r(a)(S(<T,t)+I(a,t))da, 



(2.7) 



where we implicitly assume (compare with (I2.2[) 1 that both susceptible and 
infective individuals are equally active in the population. In fact, in (|2 . 6[) we 
are also assuming that newborns are all susceptible and equally produced by 
both kind of individuals. We note that the variable 



N(a,t) = S(a,t) + I(a,t) 



satisfies the equation 



dN dN 

— + — + n{a)N{a, t) + al(a, t)=0, 
dt da 

showing how the disease-induced mortality affects the population growth. 

For the parameters regulating the infection mechanism we assume that a > 
and consider the general separable inter-cohort form of the force of infection (see 

my- 

X(a,t)=K(a) q{a)I(a,t)da , 



where the age-specific infectiousness q(-) and the age-specific contagion rate 
K(-) satisfy the following conditions 

g(-), AX)ei°°(0,at), 

q(a), K(a) > in [0,a f ] . 
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3 Well-posedness 



Existence and uniqueness of a solution to problem (12.61) can be proved following 
the standard techniques presented for example in [13] (see also [17]). For the 
sake of completeness here we give a sketch of the proof without going into 
technical details. First we denote 



B(t) = S(0,t) , 
so that, integrating S along the characteristics, we obtain 

_ f S (a - t)e~ /o [A(a-*+r,T)+M(a-t+r)] dr foj . ffl > ^ 



(3.1) 



5(0, *) ^ B ^ _ a ^ a y- f» X(a,t-a+a) da foj . ; > ft 

Moreover, using (|3.2|) and integrating I along the characteristics, we have 



(3.2) 



I{a,t) = 

e -Qt-/ V(a-t+r)rfr _ ^ 



+S'o(a - t) /J A(a - t + r, r)e QT -/o r <i£ dr j for a > i , 

n(a)e~ aa B(t - a) J Q a A(cr, t - a + o)e aa -$o <^ for i > a . 



(3.3) 



W(t) = / q(a)I(a,t)da, 



Then, if we denote 

we obtain that B(t) satisfies the following integral equation (see (|2.6jl ): 
B(t) =R d <P(Q(t)) j jT P(a)7r(a) (V aa 

+a P e' a{a -^-io K( P )w(t-a+ P ) dp da ^ B{t _ a) da 



(3.4) 



'0 



S (a - t) ( e -"t-/o M(«-t+T) da 



(3.5) 



+ ae~ ^( a -*+ CT ) da / e -"(.t-<r)-fo K(a-t+r) W(t) dr d(J 



da > , 



+ I (a- t)e- at -ti 



where all the functions are extended by zero outside the interval [0, aft\ and Q is 
defined in (|2.7p . To solve our problem, we first consider Q and as two given 
continuous, nonnegative functions. It is well known that equation (|3.5[) . being 
a linear integral equation of Volterra type with a nonnegative kernel, admits 
a unique continuous and nonnegative solution. We denote it by B(t, Q,W) to 
show the dependence on Q and W. 
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Now, for a fixed T > 0, we consider the space E = C ([0,T]; (L^O,^)) 2 ) 
and the closed set 

K ={k e E | fc(f) = (s(a, t), i(a, t)), s(a, t) > 0, i(a, t) > 0, 
ll*(-,*)IUi(o,at) < Af, ||i(-,i)|| L i(o,„t) < M} , 

with 

M = (1 + aa t )e <ll0ll,o O(o , at) (iW)T ||5o + /o || il(0 at) . 
Then, given k G S, fc(f) = (s(a, i), i(a, i)), we set 



= / r(a) (s(a,t) + i(a,t)) da , 
Jo 



W(t) 



(7(a) z(a, i) da , 



s(a,t) 



i(a,t) = < 



and define the map 

where T(k) = (s(a,t),i(a,t)) and 

fS' ( a _t) e -/o[^( a - t + CT ) w ( CT )+Ma-t+<x)]d ( x for a > t, 
[B(t - a, Q, W)ir(a)e- So K(a)W(t-a+a) da {oTt > a ^ 
' e -at-f*»(a-t+a)dv (/ ( a _ f) + S (a - t)- 

• /„* if(a - f + cr)VK(cr)e Q,T --''o < ' K"(«-t+T)W(r) dr^ for a > ^ 
7r(a)e- Qa B(t-a,Q,VF)- 

k • J Q a Jf(tr) W(* - a + ^e^So K(r)W(t-a+r) dr d(J ioT t > a _ 

Standard estimates allow to prove that T sends /C into itself and that a 
suitable power of T is a contraction, so that there exists a unique fixed point. 
Then, we can state the following result: 

Theorem 3.1. Let(So,Io) £ (^(0, a))) 2 , then there is one and only one k 6 /C , 
k{t) = (S(a,t),I(a,t)) such that 

{So(a - t) e -Sl[K(a-t+a)W( a )+Ka-t+a)\d a j Qr Q > ^ 

[B(t - a, Q, W)ir{a)e- So K(a)w(t-a+a) do f ort >a, 



S(a,t) 



(l (a-t) +S (a-t)- 



^ — at — f^ fj,(a—t-\-a) da 

/(a, t) = < ■ Jo K (a~t + o-)W(a)e aa So K(a-t+r)w(r) dr d(J ^ f ora >t, 
Tr(a)e- aa B(t-a,Q,W)- 

k • J o a K(a)W(t - a + o)e a °-So K(r)W(t-a+r) dr d(J j or t > a ; 



r 

Q(t)= r(a)(S(a,t)+I(a,t))da, 
Jo 



W(t) = / q(a)I(a,t)da. 
Jo 

Moreover, the following properties hold: 
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(i) Iim h _y ^[S(a + h,t + h) - S(a,t)] = -[K(a)W(t) + n(a)]S(a,t) a.e. in 
[0,ot] x M+ ; 

(ii) i [J(a + h, t + h) - I (a, *)] = X"(o) W(t)5(o, *) - [/i(a) + a] J(o, i) 

a.e. in [0, a f ] x M+, 

r<«; l|5(-,«)|Ui(o,-t) < (l + aa t )e fl ° d|,/3|l ^'o.„t ) (iW)T ||5o + /o|Ui(0at); 

H l|J"(-,*)IUvo,«n ^ (l + «o t )e flSll/3|li -(°- t ' (lW)r ||5o + /o|U 1 (o,«t), 
(vj there exist two constants C\,Ci > depending on M and T such that 

||*(t) - fcCtJU^co.at))* < CTie^lK-Sb./o) - (^0,-?0)||<ii(0,at)) a 
where k(t) is the solution relative to the initial datum (Sq,Iq). 



4 Search for endemic equilibria 

We now consider the problem of existence of steady states for the system (|2.6[) 
Consequently, we are concerned with the problem 

_ = -(A»+ M ( a ))S», 
dl* 

— = A* (a) S*(a) - ( M (a) + a) I*(a) , 

S*(Q) = i?o$(g*) / (3{a) (S*(a) + /* (a)) da , 
Jo 

U*(o) = o, 



where 



A* (a) = if (a) / q(o-)I*(o-)da, 
Jo 

Q* = / r(a)(S*(a)+I*(a))da. 



It is easy to see that 



S* dfe (a)=N*(a), r dic (a)^0 



(4.1) 



(4.2) 



is the disease-free equilibrium, provided that AT* (a) is given by (|2.4I) . Then we 
concentrate on the search of endemic states, that is nonnegative solutions for 
which I* (a) does not vanish identically. Our aim in this Section is to reduce the 
solution of this problem to the solution of a system of equations of the scalar 
variables 



B* = 5*(0), 
W* = 



q{a) I* (a) da. 



(4.3) 
(4.4) 
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To simplify, we denote 

L{a) = / K(a)da (4.5) 



and, integrating, we obtain 

S*(a) = B*e- WL(a K{a), 



I*(a) = B*W*n{a) I K(a)e~ w " L( -' T ^^ a -^ da . !4 (> ' 



We define the following functions: 

F{a, W) = p{a)Ti{a)^e- aa + a e- WL ^- a(a - p U(^dn , 

G{a,W) = r(a)K(a)[e- aa + u e- WL{p) - a[a - p) df^jdo , 

H{a,W) = / q{a)n{a) / K(p)e- WL ^- a( -"- p Upda . 
Jo Jo 

We note that, integrating by parts, F(ce, W) and G(a, W) can be written as 
F(a, W) = J* l3(a)ir(a) (e~ WL ^ + W J° K(p)e- WL ^- a ^- p Up^j da , 

G{a, W ) = f Q r(a)w(a) ^e~ WL[a) + W J K{p)e- WL{p) - a{ -°- p) dp\ da . 

(4.7) 

Substituting (j4~6]) in (l4~4|) we get 



,t 



W*=B*W* q(a)n(a) K(a)e' w " ' H")-^-") da da tAoS 
Jo Jo I 4 - 8 ) 

= B*W*H(a,W*) . 

Substituting (|4.6|) in (|4.1|) . we have 

Q* =B* r(o)7r(o) ^ e - w * L(a) + W* K(a)e- w ' L(a) - aia ~ a) da^ da 
and, using (|4TT)) and fO]) . 

H if(a,W*)- 1 j 

Now our problem has been reduced to solving the following system of equations: 

j^* _ 

H(a,W*)' 

/G(MnN ^ 
\H(a,W*)J v ' y 



The second equation in (|4.10[) depends only on W* : once we have a solution of 
it we can substitute it in the first equation. So, in order to simplify the second 
equation, we write 

<p{a, W) = Rfa WO (4-H) 

and we are left with equation 

<fi(a,W) = 1. 

We have that 



.it 

^(a,0) =i?^$ ' ' 



J" r(a)ir(cr) d a 
H(a,0) 



and 



lim <p(a, x) = , 



so the second equation in (|4.10|) has at least one solution W* if y>(a, 0) > 1. 
Since $(•) is decreasing in [0, this latter condition is equivalent to 

Jo at r(a)7r(g) rig 
iT(a,0) 

where is the solution of (|2.3[) , that is to (sec (12.41) ) 

J&gfa °) = N *(p)H{a, 0) > 1 ■ (4.12) 
J a r(o>(o-) da 

Condition (14. 12|) can be written as 

fl§>l, (4.13) 

where i?§ is the basic epidemic reproduction ratio, i.e. the number of secondary 
cases which one case would produce in a completely susceptible population. This 
number is defined in this context as the spectral radius of the next-generation 
operator Q (see [8], Chapter 7) 

(Qu)(a) = / g(a,a)u(a) da , 
Jo 

where 

g(a,a)= f K(a)q(p)N*(a)^\e- a ^ dp . 

Jo- ^(a) 

It is easy to see that 

m= [ K(a)N*{a) f g( CT )Z!M e -°(»-<0 da da . 

J0 Ja ""(a) 
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Since N*(0) = ^.ff for any a S [0, at], we can write 



,.t 



N*(0)H(a,0) = [ q(a)n(a) [ K(a) e - a ^- a) ^$dada 

7r(a) 



o Jo 



K(a)—^r / q{a)-K{(j)e- a[ ' 7 - a) da da = R L Q . 

1"W Ja 



Thus (|4.13p is a sufficient condition for the existence of an endemic equilib- 
rium. 

Next section is devoted to discuss uniqueness of equilibria under such con- 
dition. 

5 Discussing uniqueness of an endemic equili- 
brium 

Condition (|4.13p is sufficient to have at least one endemic state, but for unique- 
ness we need some additional assumptions. A first uniqueness case occurs when 
the disease does not induce mortality (see also [6]). 

Theorem 5.1. If a = and Rq > 1, there exists one and only one endemic 
equilibrium of the problem 



Proof. Existence of endemic equilibria is equivalent to existence of solutions of 
the system (|4.10|) . Since the function 

r(cr)7r(cr) da 



j; q (aMa)f a K(p)e-WL( P ) dpd(J 

is increasing in W and tends to infinity as W tends to infinity, there exists W 
such that 

I if r(a)ir(a)da \ f a * , ^ , , 

<p(0,W) = $ —i ^ / p( a )Tv{a)da 

V Jo 9^Ma) J° K{p)e-wH P )dpdo j h 

is decreasing for W in [0, W], with ip(0, W) < 1, and non-increasing for W > W. 
Now, Rq > 1 implies <p(Q, 0) > 1 and so there exists one and only one solution 
oftp(0,W) = l. □ 

The proof of another condition for a > 0, requires the following Lemma (see 

El) 

Lemma 5.2 (Hadeler and Dietz). Let g, u and v be locally summable real 

functions on the (finite or infinite) interval (a,b), with u(x), v(x) > 0, such 
that 

x < y =>■ g{x) < g(y) and v(x)u(y) > u(x)v(y) . (5-1) 
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Then 

r-h r-b f-b f-b 

I g(x)u(x) dx ■ / v(x)dx> / g{x)v(x)dx- / u(x) dx , (5-2) 

J a J a J a J a 

provided the integrals exist. 

Moreover, when u, v and g are positive and is strictly increasing, then 
i5.2\) holds with strict inequality. 



Then we obtain the following result. 
Theorem 5.3. If for any < p\ < pi < a* , 

q(a)n(a)e- a,T da / r(a)n(a)e~ a,T da 

P2 ♦ (5-3) 

<K(p 2 ) r(a)n(a)e- " 7 da q(a)Tt{a)e- aa da . 

J pi J p2 

and Rq > 1, then there exists one and only one endemic equilibrium of the 
problem h2.b}) . 

Proof. Fixed a > 0, Rq > 1 implies (f(a, 0) > 1. By definition F(a,W) is 
decreasing in W. Then, if we prove that ^"'j^j is increasing in W, since 
\imw^>+Qo g^vi?) = +oo, we have that ip(a, •) is decreasing on a suitable inter- 
val [0, W], where <f> ( G ( a ')P ) > o, and this gives our claim. Since 

_d_ ( G(a,W) \ = d J^lH{a,W)-G{a,W) d -^l 
dW\H(a,W)J (H{a,W)) 2 

we study the sign of the numerator GyyH — GHw- If we denote 



and 



u{p) = K(p)e~ WL{p)+ap / q(a)Tr(a)e~ aa da 

V( P ) = e ~WL(p)+ ap / r((7)7r((j)e - QCTd(T ^ 



we get: 



G W H — GH W = - a / L(p)v(p)dp I u(p) dp 
Jo Jo 

+ a v(p) dp / L(p)u(p) dp 



+ / r{a)ir{a)e- aa da / L(p)u{p) dp 
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Now we apply Hadeler-Dietz's lemma. Since L(p) by definition (see (|4.5[) ) is 
increasing, we choose g = L. The second part of condition (|5.ip coincides with 
(|5.3I) . so we obtain that 

G W H -GH W > \ r(a)Tr(a)e- a,T da / L(p)u(p) dp > 0. 
■/o Jo 

□ 

Remark 5.4. Condition i5.3\) is verified for example when there exists a func- 
tion h(-) such that q(a) = h(a)r(a) and both K(-) and h(-) are non- decreasing 
functions. 

We analyzed so far sufficient conditions to have uniqueness. We can also 
show examples in which multiple equilibria occur. The following case corre- 
sponds to a very special situation but it helps to understand the mechanisms 
responsible for non-uniqueness. 

Let us consider the following choices 

o+ = f r(a) = 1 

a = 10 p(a) — tana 

R$ = 27 K , , f 1 ifae[0,|]u[|,|] (5.4) 



[3(a) = 1 K{a) \ ifae (5,f) 

g(a) = 1 = max{l - ^,0} 

Then wc have 



6 

7f(a) = cos(a) , L(a) = ■( | ifae(|,| 

a _| ifoG[f,f 

Thus we may explicitly compute the functions F, G and -ff and we find that 
the function (see (|4.11|l ) 



F(10, W) 



has the graph shown in Figure [T] This leads to the following result 
Theorem 5.5. Problem \2. 6|) with the choices j5.4\ ) has two endemic equilibria. 
As an extension of this special case we consider the following forms 

r q = ^Y and ®( x ) =max|l- ^-,oj (5.5) 

and in Figure [2] we show the endemic equilibria curves versus the parameter a, 
for different choices of X. The special case of (|5.4[) corresponds to X = 18. The 
behaviour becomes more evident as X increases. 
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2 4 6 8 

W 



Figure 1: the graph of (p(10, W) for the choices (|5.4[) . 




a 



Figure 2: endemic equilibria curves for X = 10,14,18,22,26,30,34,38,42 in 
<|53]) . 



6 Characteristic equation 

In order to investigate the local asymptotic stability of the equilibria we consider 
the characteristic equation. For this standard technique see for example [13] . 
In the following sections we assume that $ is differentiable at Q* d (in the case 
of the disease free equilibrium) and in Q* dehned in (|4.ip (in the case of the 
endemic equilibria). 

A convenient way to represent problem (|2.6|) in an equivalent form is to use 
the variables B{t) as in (|3~Tj) and W(t) as in ([3^4]) . In fact, using (|3~2|) and ([3~3]) . 
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we get a system of two integral equations that, for t > , is given by 

B(t) = ^ jT r(o)7T(o) (e~ R K(a)W(t^a+a) da 

+e- aa J^K(<j)e°" T -J'» K(p)w{t-a+ P ) dp W ( t _ a + a ) dajB(t - a) da)j- 
•jf /3(a)vr(a) (e~ *C")W(t-«+«0 da 

+e- aa J^K(a)e a,y '^ K(p)w(t-a+ P ) dp w( t _ a + CT ) d(j^jB(t - a) da , 

W(t) = / q{a)ir{a)e- aa K(a)e aa -fo K{ P )W(t-a+ P ) dp w ^ _ a + CT ) da . 
Jo Jo 
■ B{t — a) da . 

(6.1) 

Actually, the limiting equations of the two general integral equations for B(t), 
W{t) coincide with (I6.1|) . Moreover, the constant solutions Bit) = B* , W(t) = 
W* of (|6.ip correspond to the different steady states of the problem. Namely 
we have: 

• the trivial solution B* = W* = 0; 

• the disease-free equilibrium corresponding to W* — and B* provided by 
the equation 



1 = i^B* r{a)ir(a) daj ; 



(6.2) 



the endemic equilibria provided by system (|4.10p . 



Note that the case of the disease-free equilibrium corresponds to (|4.2j) with (|2.4I) 
and (compare (|6.2p with 



(6.3) 



Q* d = B* I r(a)Tr(a)da. 
Jo 



Finally, the endemic equilibrium is given by (|4.6j) . 
In order to linearise (|6.1j) at (B*,W*), we let 

b(t) = B{t) - B* and w(t) = Wit) - W* . 

Then we obtain the linear system 

b(t) = / *x(a)6(t — a)da+ / ^ 2 {a)w{t - a) da , 
Jo Jo 

wit) = / * 3 (a)6(t-a)rfa+ / * 4 (a)u;(i - a) da , 
Jo Jo 



(6.4) 
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where the convolution kernels ^i, ^2, ^3 and ^4 are given by 

■ tt(o) ^e-^ L W + jT A-(a) e -"( Q - CT )- w * L ^)^ , 
* 2 (o) = - aB* £ Ufa(Q*)P(a) + B*^^r(a^j 7r(a)K(a a)- 



e -<*(°-p)-W'L(p) dpdat (6.5) 

¥ 3 (a) =W*g(a)7r(a) / K(a) e -^ a - a) - w " L ^da , 
Jo 

* 4 (a) =£* / g(cr)7r(cr)/s:( < T - a)- 

J a 

. ^-aa-iy*L( CT -a) — V^* J K ^e^- p) ~ W * L(p) df^j da , 



with Q* defined in (|4~Ti 

Taking Laplace transforms in (16. 4[) . we obtain the following characteristic 
equation (see [13]): 

¥(A) = (1 - $i(A))(l - $ 4 (A)) - $ 2 (A)$ 3 (A) - , (6.6) 
where for i = 1, 2, 3, 4 

/•+0O 

*i(A)= / e- AQ ^(a)rfa 
Jo 

represents the Laplace transform of ^ (a) . 

To study the stability of an equilibrium, we need to determine the location 
of the roots of (|6.6|) . In the following section we will be concerned with the 
disease-free equilibrium, then we will investigate the role of the parameter a. 

7 Stability of the disease-free equilibrium 

In the case of the disease- free equilibrium (|4.2p , the convolution kernels become 
*i(a) =^(o)tt(o) + <B*$'(Q*) r(o)7r(o) , 

* 2 (o) =B*(e~ aa - 1) / P(<t)tt((t) K(a - a) da 

J a 

+ R d {B*) 2 &{Q* d ) / /3(a)Tr(a) da{e~ aa - 1) r{a)n{a)K{a - a) da , 

JO J a 

*s(a) =0, 

* 4 (a) = B*e~ QQ / g(cr)7r(CT)i4T(cr - a) da , 
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where Q* d is given in (|6.3p . 

Hence, the roots of f|6 . 6[) are the union of the roots of the two equations 

$i(A) = 1, * 4 (A) = 1, 

that can be considered separately. Note that the first of this equations is ex- 
actly the characteristic equation (|2.5p for the demographic problem (|2.ip in the 
absence of the disease. Thus we have 

Theorem 7.1. If Rq > 1, the disease-free equilibrium is unstable. If Rq < 
1, the disease-free equilibrium is stable or unstable depending on whether the 
equilibrium \2.J$ of the population is stable or not for (2.1\) . 



Proof. Let R% > 1. We have that 

pa* no^ par 

$4(0)=/ $4{a)da = B* / e~ Qa / q(V)7r(cr)Jf (cr - a) da da 

Jo Jo J a 

= B*H(a,Q) . 

Since 

G(a,0) = Q* d n * 

H(a,0) N*{0)H(a,0) Vd 

and <&(■) is decreasing on [0, Qj], one has 



Now, comparing (|6.2p with the previous equation, we get that 

1 



B* > 



H(a,0) 



and consequently $4(0) > 1. This means that 1 — i&4,(t) has exactly one root 
on the positive real line. 

If Rq < 1, then $4(0) < 1 and, if A G C has positive real part, 

1 - $ 4 (A)| > 1 - 



/ e~ Xa ^i(a)da 
Jo 

>1-/ *4(a)da = 1 - * 4 (0) > 0. 
Jo 



So, there are no zeros of 1 — ^(A) with positive real part. The study of the zeros 
of 1 — \l/ 1 (A) coincides with the stability analysis of the non trivial equilibrium 
of the total population (equation (|2.ip ) and this concludes the proof. □ 

From the previous Theorem we see that the epidemic reproduction ratio i?g 
determines also stability of the disease free equilibrium. In the next Section we 
discuss how instability may depend on the disease induced mortality a. 
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8 Destabilizing effect of the extra-mortality 

Since an analytic study of the characteristic equation at the endemic equilibrium 
does not seem to be possible in the general case, we present a specific example. 
We consider a special set of parameters that shows how the extra-mortality a 
can destabilize the endemic equilibrium. 

Let us first consider the following choices: 

7T 3 3 

a) = — , (3(a) — — sin(2a) , r(a) = — sin(2a) , p(a) — tana . (8.1) 



With these choices and a = the convolution kernels (|6.5[) for the endemic 
equilibrium become 

*i(o) = (1 + R$B*&(Q*)) | sin(2a) cos(a) , 



*2(a) =0, 

* 3 (a) =q(a)co S (a)(l-e- w * L ^) , 
,«+ 

* 4 (a) / q(a) cos(a-)^(c7 - a)e" w l(ct W , 

J a 

where B* and Q* are given by (|4. 10[) and (|4.9I) . Then we have the following 
result. 



Proposition 8.1. Let \8.1\) be satisfied. If a = cmd 

fj&B*$'(Q*) = -5, (8.2) 

then the characteristic equation H6.6}) has two imaginary roots X± = ±5i and 
any other root has negative real part. 

Proof. Since ^2 (a) = 0, the characteristic equation is simply 

*(A) = (1-$!(A))(1-$ 4 (A)) = 0. 
If A has nonnegative real part, we have 



R(*4(A)) 



,t „„t 



B* e- Xa q(a)n{a)K{a - a)e~ w ' L( - a) da da 

J0 J a 

< B* / q(a)n(a)K(a - a)e~ w ' L ^ da da . 

J0 J a 



Since 



f f q(a)ir(a)K(a - a)e~ w * L ^ da da = f q(a)7r(a)e' w " L{a) f K{p)dpda 

Ja Jo Jo 

< r q {a)^{a) ( a K(p)e- w * L M dpda = H(0,W*) = , 
Jo Jo u 
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li 1 



< 1. Then the roots of the characteristic equation with nonnegative 



(*4(A)) 

real part can be found only if \&i(A) = 1, i.e. 

(! + T )o / e_AQ sin(2a) cos(a) da = 1 (8.3) 

with t — RqB*&' (Q*). We look for values of r < for which a couple of imag- 
inary roots A = ±iuj of (|8.3[) exist. This is equivalent to solving the following 
system: 

{Jq sin(a;a) sin(2a) cos(a) da = , 
(1 + r)| cos(wa) sin(2a) cos(a) da = 1 . 

The first equation gives 

sin — L r 5 — - sin — ^-s — 1 sin — ^ — - sin — ^ — 1 4w cos ^ 

2 i 2 2 2 2 n 



4(w-3) 4(w-l) 4(w + l) 4(cj + 3) (cj 2 - 9)(w 2 - 1) 

and the solutions are = ±(2fc + 1) with = 2, 3, Correspondingly, the 

values of r = 77. are given by the following equality 

2 

T fe = — 1 

3 J Q 2 cos((2fc + l)a) sin(2a) cos(a) da 

= - ? ^-1, 

cos(2(fc-l)a) , cos(2fca) _ cos(2(fc+l)a) _ cos(2(fc+2)a) 
fc— 1 ' fc fc+1 fc+2 

which means 

I 1-4fc for fc even , 

Tfc = ' 



-4fc -8fc-3 



for fc odd . 



I 3 

It is easy to verify that is decreasing with fc. Then T2 = — 5 is the maximum of 
the Tfc's. Since by the implicit function theorem one obtains that ^r(-5) < 0, 
A = ujfi = ±5i are the first roots to cross the imaginary axis and this completes 
the proof. □ 

We now let a become positive, and analyze how the stability of the endemic 
equilibrium changes with a. For simplicity we make a further choice considering 

$(s) = maxjl — — ,0| , q(a) = - sin(2a) , K(a) = a. (8.4) 

By ()8.1j) and Remark l5.4l for every fixed a there exists one and only one endemic 
equilibrium. So, equation ip(a, W) — 1 has only one solution (0, Wq), with 
a = 0. With the choices made in (jgTTJ) . and (|K4l . we get 

i? d = 6 and ff(0,W o *) = ^. 
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so that Wq ~ 5.04512. Moreover at this point we have °- Thcn 

for a > sufficiently small, there exists a function W*(a) such that 

^(a,VF*(a)) = 1 (8.5) 

and we obtain a branch of endemic states (B*(a), W*(a)). 

Now we show that this endemic equilibrium changes its stability, in the 
sense that at a = the roots of the characteristic equation cross forward the 
imaginary axis. 

We write A = £ + iuj, so that we can consider the characteristic equation 
(|6.6I) at the equilibrium (B*(a), W*(a)), as a system depending on a 

Fi(a,C,w) = m(X) = 0, 
F 2 (a,C,w) =9f*(A) =0. 

By Proposition O we know that Fi(0,0,±5) = F 2 (0,0,±5) = 0, while 
detailed calculations show 

^(0,0,±5) = --^(0,0,±5), 

fit? rip 
^(0,0,±5) = ^(0,0,±5), 

so that the jacobian of the system with respect to £ and uj at (0, 0, ±5) is equal 

to 

2 / N 2 



v (?c -(0,0,±5)j +^(0,0, ±5; 

and by numerical computation (the evaluations are obtained using Mathemat- 
ica, www.wolfram.com) we conclude that it is positive. Then for a > suf- 
ficiently small there exist ((a) and ui(a) such that £(0) = 0, w(0) = ±5 and 
Fi{a, (,{a), 00(a)) = for i = 1, 2. Concerning the sign of C'(0), again by numer- 
ical computation we have 

C , (Q) = = gjO; 0, ±5)ffi(0, 0, ±5) - ffi(Q, 0, ±5)^(0, 0, ±5) > q 
(^(0,0,±5)) 2 +(^(0,0,±5)) 2 

and we conclude that, in this example, disease-induced mortality yields insta- 
bility. 

The calculations to obtain the conclusion above are standard, but for the 
reader's convenience we give below some details. With the choices made in 
this special case we have the following kernels where we have highlighted the 
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dependence on a 



= ( 6<&(Q») - 1Q Q(Q°l a)) ) |sin(2o)cos(o)- 



o - a (a-a) ±—o d(J ^ 



* 2 (a,o) = -aB*(o)(6$(Q*(a))- 



im(Q*{a))J 2 

W(a) „2 



sin(2cr)cos(a)(a-a) / e ~ a ^-P) ~P dpda , 



a — a 

" a) -2 



$ 3 (a,a) =Ty*(a)^sin(2a)cos(a) / ae-^ "^- 1 ^^^ 

2 Jo 

3 /"a 

^(a, a) =S*(a)- / sin(2cr) cos(er)((7 — a)- , 

2 Ja 

W*(c«) 2 f a , ^ W*(a) 2 

e = CT -a / e - Q ( ff -") r^'dp da. 



a - — a 



Then we evaluate the following derivatives (they are the only derivatives we 
need) 

-^(0, a) = sin(2a) cos(a) l^-3W*(0) + 6a - 6 y e —° da J , 

p r w(o) 2 

— — (0, a) = 50 / sin(2a) cos(cr)(cr - a) / e 2 p dpda, 
dot J a J a _ a 

where we have used 

W(o) = (™-w(o)-i) gg^W 
1 J Ub ( j 7§(o,r(o))' 



obtained by Dini's Theorem applied to (|8.5I) . The above expressions are used 
to compute the following quantities 

A = 1 - 2 cos(5a)* 4 (0, a) da ~ 0.83432, 
Jo 

73= / 2 sin(5a)* 4 (0,a)da ~ 0.236031, 
Jo 

C = / cos(5a)— -!-(0,a)da ~ 2.52632, 
Jo oa 

73=/ sin(5a)— ^-(0,a)da ~ -0.397545, 
Jo oa 

E= I 2 cos(5a)* 3 (0,a)da ~ -0.147547, 
Jo 

F = 2 sin(5a)* 3 (0,a)da ~ -0.177344, 
Jo 
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G = 



I 



■2 



cos(5a) 



d^ 2 



(0,a)da ~ -0.0744986, 



da 




9*2 



(0,a)da ~ -0.249824, 



Jo 

that finally allow us to evaluate 



da 



-^(0,0, ±5) 




AC - BD - EG + FH ~ -1.98062, 



± (-BC + AD + FG + EH) ~ ^0.877897, 
|U+|b~ 0.5669, 

± ( t^-B - §A ] ~ tO.440352, 




so that C'(0) ~ 1.42878 > 0. 

9 Numerical exploration 

In the previous Section we have produced an example showing that the parame- 
ter a, representing disease-induced mortality, can actually modify the dynamics 
of the system and that periodic solutions are possible via Hopf bifurcation. In 
order to explore the model in a systematic way, we now resort to a numerical 
method that allows to determine the roots of the characteristic equation and 
follow their displacement as a varies. The method, proposed in [3], provides 
numerical approximations to the rightmost part of the characteristic spectrum 
associated to the model linearized around the equilibrium to be investigated. It 
is indeed well-known that the zero solution of this latter is asymptotically stable 
if and only if all the characteristic roots have strictly negative real part. 

The numerical scheme developed in [3] is actually devoted to the stabil- 
ity analysis of the scalar Gurtin-MacCamy model [10], but it can be extended 
straightforwardly to the m-dimensional system (to > 1) 



where P : [0,a^] x [0,+oc) — > W n is the m-vector of population densities, 
M,B : [0,a^] x R" — > M mxm arc the matrices of mortality and fertility rates, 
respectively, and S : [0, +oc) — > R" (n > 1) is the n-vector of population sizes, 



( dP dP 



< 




(9.1) 



LP(o,0) = P o (o), 
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i.e. a selection of n homogeneous population sub-classes through the weight 
function Q : [0,at] -> M" xm . The epidemic model fl22J} we are interested in fits 




Figure 3: bifurcation diagram of equilibrium W* as a varies (X = 34 in (|5.5 



into (|9.ip by choosing m — 2, n = 2, the population vector 



P(a,t) = 0S(a, 



and the size vector 



s(t) = m),w(t)y , 

which corresponds to selecting 



Q{a) = 



r(a) r(a) 
q(a) 



Moreover, the matrices relative to the vital rates are given by 



M(a,S(t)) = 



K(a)W{t)+n(a) 







-K(a)W(t) a + /j,(a) 



and 



B(a,S(t)) 



R$P(a)*(Q(t)) R d Q P(a)<S>(Q(t)) 











Within this framework we consider the special choices of Section [5j namely 
(|5.4p and (|5.5[) with X = 34, letting a vary along the corresponding equilibrium 
curve represented in Figure [5] (the third curve from the right). For each point on 
the curve (included those corresponding to the disease-free equilibrium for W* — 
0) we are able to compute the rightmost characteristic root (rounded to machine 
precision, see [3]) and thus to say whether the corresponding equilibrium is 
locally asymptotically stable or not. The overall situation is illustrated in Figure 
[3] where solid lines denote stability and dashed lines denote instability. 
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We start our analysis from the right hand branch of the disease free equi- 
librium by investigating, for instance, the spectrum for a — 24 (A in Figure [3] 
and Figure 01 : the rightmost roots have negative real part and hence the trivial 
equilibrium is stable. By decreasing a the real root moves to the right until at 
a — ai ~ 20.1143 it crosses the imaginary axis rightward (B). This first bifur- 
cation makes the trivial equilibrium lose its stability and a branch of endemic 
equilibria raises. Following the disease-free branch by further decreasing a it is 
confirmed that the instability persists (C, a = 10). 

Going back to the bifurcation point at ai, we now follow the endemic branch 
by increasing a. The branch is unstable (D, a = 21) until at a = «2 — 
22.8495 the leading root crosses the imaginary axis leftward (E). This second 
bifurcation makes the endemic equilibrium gain back its stability. Then the 
branch continues by decreasing a again and at a = 21 (F) the second endemic 
equilibrium is stable opposite to the first one for the same value of a (D). The 
branch remains stable by further decreasing a (F, a = 10). 

Similar trends are obtained for other values of X for which double endemic 
equilibria exist. On the other side, for those values of X for which only one 
endemic equilibrium exists (for instance X = 10, 14 in Figure [2]), it can be 
observed by the roots computation that at the first bifurcation the disease-free 
equilibrium loses its stability in favour of the endemic one. This latter then 
preserves its stability for decreasing values of a down to a = 0. 

As a consequence of the destabilisation shown in Section[Sl it is possible that 
under different choices of the parameters a Hopf bifurcation occurs. In fact, if 
we consider the same choices made in (15.41) for and K, but set 



/3(a) = | sin (2a), u (a)-taiia 

q(a) = sin (2a), */• { r\ * m ( 9 - 2 ) 

) \ 3 ■ rn \ $(x) =max{l- ^,0}, v ; 
r(a) = § sin (2a) , w 18 5 J 5 

we get the bifurcation diagram represented in Figured] (left). The black dot, 
corresponding to a ~ 0.6743, indicates a Hopf bifurcation through which the 
equilibrium on the stable branch (as usual in solid line) loses its stability and 
a limit cycle arises. The right figure shows the existence of the corresponding 
couple of characteristic roots crossing the imaginary axis from left to right as 
a decreases. Thus, in this case the disease-induced mortality has a stabilising 
effect. 

We can also show a choice of the parameters for which the diagram is closed 
and there is, for increasing a, a destabilisation and, successively, a stabilisation. 
In fact, for 

3 

q(a) = Wa, r(a) = - sin(2a) , X = 15 and Jig =1.35, (9.3) 
o 

one has the bifurcation diagram of Figure |6] 
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Figure 4: relevant characteristic roots for different equilibrium points, referring 
to the bifurcation diagram of Figure [3] (X = 34 in (|5.5|l ). 
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Figure 5: bifurcation diagram of equilibrium W* as a varies (left) and relevant 
characteristic roots (right) for a ~ 0.6743, value at which a Hopf bifurcation 
occurs (black dot in the left figure). 



5 0.5- 



Figure 6: bifurcation diagram of equilibrium W* as a varies. The dotted line 
corresponds to periodic solutions. 



10 Conclusions 

The epidemic model studied here differs from the classical ones for the pop- 
ulation size, which is not constant. As pointed out in [12], this leads to the 
possibility of having multiple endemic equilibria. In our case, which is reg- 
ulated by a very different mechanism, we are able to show that under some 
conditions two endemic equilibria occur. One of them is unstable and the other 
is stable. 

It is well known that, when the population size is fixed and transmission is 
inter-cohort, one has uniqueness of endemic states. See for example [I] and [5]. 
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The stability of the steady states is studied in [IS], [7] and we obtain here an 
analogous stability change for the trivial equilibrium. 

For the sake of simplicity we have limited ourselves to investigating the 
case of inter-cohort transmission. Allowing for a general transmission kernel 

\(a,t) — Jq K (a, a) I (a, t) da makes the analysis of the model very difficult 
already in the case of constant population size (see |14j). while we wished to 
emphasise the peculiarities due to infection-related deaths. 

The situation where multiple endemic equilibria exist may not be realistic. 
In fact, the model studied assumes that contagion can happen only in the first 
and in the last period of an individual life. A somehow similar situation was 
studied in [5], where the infection rate is piecewise constant. More precisely, 
juveniles and noncore adults cannot be infected, while only core adults can. 

We have shown that the model with age-structured contact rates and variable 
population, because of infection-related deaths, has a much richer bifurcation 
diagram than models with only one of these features. For example, if (|8.ip , 
(|8.2I) and (|8.4[) hold, there is a change of stability and the presence of the extra- 
mortality destabilises the endemic equilibrium. Moreover, if either (I9.2[) or (|9.3p 
holds, Hopf bifurcations occur. 

We believe to have shown the main theoretical steps necessary to determine 
the stationary solutions of the model and their stability properties, adding the 
use of the numerical tool provided in 3] whenever analytical conclusions are 
hard, when not impossible, to draw. We also have highlighted the complexity 
of this model, that presents different kinds of dynamics. 
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